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Abstract 

We consider a vertex model in d dimensions characterized by lines which 
run in a preferred direction. We show that this vertex model is soluble if the 
weights of vertices with intersecting lines are given by a free-fermion condition, 
and that a fugacity —1 is associated to each loop of lines. The solution is 
obtained by mapping the model into a dimer problem and by evaluating a 
Pfaffian. We also determine the critical point and the singular behavior of 
the free energy. 

PACS number: 05.50.+q 
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I. INTRODUCTION 



An outstanding unsolved problem in statistical mechanics has been the exact solution of 
lattice models in three or higher dimensions. While a host of lattice models has been solved 
in two dimensions [|T|, very few exact results are known for three or higher dimensions. 
Some two decades ago Suzuki |^ solved a three-dimensional Ising model with pure 4-spin 
interactions, which turns out to be a two-dimensional model in disguise. Similarly, a three- 
dimensional 5-edge model solved by Orland is equivalent to a two-dimensional dimer 
problem The first exact solution of a genuine three-dimensional model was that obtained 
by Baxter |^ for the Zamolodchikov model |^, a 2-state interaction-round-a-cube (IRC) 
model involving some negative Boltzmann weights. Since then Baxter and Bazhanov have 
further extended the solution to the IRC model of states for general N. Very recently, 
we have solved a vertex model in arbitrary dimensionality p[. This is a vertex model 
described by lines which run in a preferred direction on the lattice and do not intersect. 
Like the three-dimensional IRC model of Refs. and 0, this model also involves negative 
Boltzmann weights. While the nonintersecting feature makes it possible to use this model 
to describe flux lines in type II superconductors the more general version of the model 
with vertices of intersecting lines remains unsolved. 

Some twenty-five years ago one of us and Fan |TO|,[TT]| introduced the term free-fermion 
model to lattice-statistical vertex models satisfying a certain free fermion condition. When 
this condition holds, the model can be regarded as one describing free fermions on a lattice 
|§| and becomes soluble. But soluble free-fermion models have been restricted mainly to two 
dimensions [0. In this paper we extend the solution of [§] to a free fermion model which 
includes vertices of intersecting lines in any dimension. This yields a new class of soluble 
free fermion models in arbitrary dimensionality. 
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II. THE MODEL 



In the vertex models considered in [Q], vertex configuration are described by lines em- 
bedded on a lattice subject to: 

i) lines do not intersect, 

ii) lines are oriented in a preferred direction, and 

iii) each loop of line possesses a fugacity —1. 

Now we relax condition (i) by allowing lines to intersect. For definiteness we consider a 
Cartesian lattice C in d dimensions. Place bonds along edges of C with the restriction that 
the number of bonds incident at each vertex is zero or even, with half of the bonds incident 
in the positive and the other half in the negative axes directions. The bonds then form lines 
which run along the main diagonal of C. Thus, we have a D- vertex model with D different 
vertex configurations, where 




where r is the number of lines intersecting at the vertex. This leads to D = 6,20, .... for 
d = 2,3, respectively. 

Number the negative axes incident at a vertex {n} = {1, 2, d} and the positive incident 
axes {p} = {d + l,d + 2, 2d}, such that {i, d + i} refer to the axis in the ith direction. 
Denote by uj^n'}{p'} the weight of a vertex with bonds (lines) incident at {n'} G {n} and 
{p'} ^ {p}- Foi' vertices with no incident bonds we take the weight Uq = 1, and for vertices 
with two incident bonds, one in the negative direction i and one in the positive direction j, 
we take the weight to be Uij. For vertices with four or more bonds the weights are given by 
the free-fermion condition 

^{n'}{p'} = ^{-if'' n ^ij 

P ie{"'}je{p'} 
= detl^ijl, (2) 

where the summation is taken over all pairings P of axes into pairs {ij}, i € {n'} and 
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j G {p'}, and Sp is the signature of P, the number of transposition needed to bring the 
indices of the pairing to the canonical ordering 

{l,d+l}{2,d + 2}-- -{(1,2(1}. (3) 

Here, it is understood that indices not contained in {n'} and {p'} are deleted in (3) in 
determining the canonical ordering. For d = 3, for example, (2) gives 

1^1245 = 1^141^25 ~ 1^151^24 
'^1256 = ^16^25 ~ ^^15^^26 

1^1246 = ^u^2e ~ ^iG^24, (4a) 
for {n'} = {12}, and other similar relations for {n'} = {13}, {23}, and 

1^123456 = ^14^25^36 + '^15'-^26'^34 + '-^16'-^24'^35 

— a;i4Co'26t^35 ^ ^^15'^24'^36 ^ ^^16^251^34, (4b) 

for {n'} = {n} = {123}. Note that the expression (2) can be written as a determinant as 
indicated. The model considered in corresponds to taking Uij = y/zizj, Zi = Zi+d so that 
the weights (4a) and (4b) vanish identically, indicating that vertices with 4 or more bonds 
are forbidden. 

The vertex weight (2) can be represented graphically by decomposing the vertex config- 
uration into continuous lines indicated by the pairings. That is, for each uJij one connects 
the negative axis i to the positive axis j. For the vertex weight (4b), for example, this leads 
to the graphs shown in Fig. 1. In this picture and assuming periodic boundary conditions, 
lines form loops wound around the lattice. Let i be the overall number of loops when each 
vertex configuration is decomposed into its canonical ordering. This permits us to define a 
partition function 

Z= E (-1)' n ^KHp'}- (5) 

vertex config. vertex 

Now, substitute (2) into (5) and observe that when all pairings are in the order of the 
canonical form we have 6p = 0, and whenever the pairing of two lines interchanges at a 
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vertex, such as changing {14, 25} to {15, 24} in the first hne of (4a), both £ and 5p change 
by 1. It follows that we can rewrite the partition function (5) as 



where the summation is taken over all possible loop configurations of £ loops, and the product 
is taken over all a;^ factors only, which are present along the loops. In the next section we 
show that, in the form of (6), the partition function Z is precisely a dimer generating 
function. 



Introduce a dimer lattice C by expanding every lattice point of C into a "city" , which 
we choose to be the one consisting of 2d points located on the 2d incident axes (n) and 
(p). The situation of d = 3 is shown in Fig. 2. Let N and 2dN be the respect numbers 
of lattice points in C and C . Connect as shown in Fig. 2 each point in {n\ to every point 
in {p}, and vice versa. There are now two kinds of edges in C : intercity edges connecting 
cities and intracity edges within a city. Starting from a given vertex configuration on C, 
we cover each intercity edge on C by a dimer (resp. leave it empty) if the corresponding 
edge on C is empty (resp. covered by a bond). The remaining uncovered points on C are 
then covered by dimers placed along intracity edges. This latter covering corresponds to 
the different pairings of the remaining points and is generally not unique. For example, the 
different coverings corresponding to 00124,5 and a;i23456 given by (4a) and (4b) together with 
the associated signs are shown in Fig. 3. In this way, we have mapped vertex configurations 
into dimer coverings. 

We now establish that there is a one-one correspondence between dimer coverings of 
and loop configurations in (6). Consider first dimer coverings on C Superimpose any given 
dimer covering Cj with a standard one, Co, in which all intercity edges are covered by dimers. 
The superposition of two dimer coverings produces a graph of transition cycles, or polygons 



E 




(6) 



loop config. 



III. A DIMER PROBLEM 
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p^ . In the present case, the transition cycles belong to one of two kinds: double dimers 



placed on intercity edges, and polygons wound around the lattice forming loops. Thus, each 
dimer covering is mapped into a loop configuration. Conversely, to each loop configuration 
in (6), there exists a unique dimer covering Ci which, when superimposed upon Co, produces 
the loop configuration in question. This completes the proof. 

It is well-known that dimer coverings are generated by a Pfaffian, provided signs can be 
fixed correctly [l^. In the usual dimer problem all terms in the dimer generating function 
have the same sign, and for this reason the generating function can be written as a Pfaffian 
only for planar lattices. In the present case, however, the partition function (6) is precisely 
a Pfaffian as we now see. 

Number the 2dN lattice points of C in the increasing order along the main diagonal of 
£, and within a city in the order as shown in Fig. 2. Direct all intercity edges of £' in the 
positive direction and all intracity edges of C as shown in Fig. 2. In this way all edges are 
directed in the positive directions. Further, the factor WoJij in (6) is generated by taking the 
dimer weight 1 for intercity edges and weights uJij,i < j, for intracity edges. To fix the sign 
associated with this product, we consider a typical term in the Pfaffian corresponding to the 
dimer configuration Cj. The sign of this term relative to the term corresponding to Cq is 
the product of the signs of the transition cycles produced by the superimposition of Ci and 
Cq |]13|. The rule is that each transition cycle carries a sign (—1)""+^, where n is the number 
of arrows pointing in a given direction when the transition cycle is traversed. For transition 
cycles consisting of double dimers we have n = 1 so that the sign is always positive. For 
transition cycles consisting of loops around the lattice we have, since all edges are directed 
in the same direction, n = even, so that the sign is always negative for each loop. Then the 
overall sign of Ci relative to Cq is (—1)^. This establishes that the partition function (6) is 
precisely a Pfaffian. 
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IV. EVALUATION OF THE PARTITION FUNCTION 



We now evaluate the PfafRan. 

Let N — Ni X N2 X ■ ■ ■ X Nd, where Ni be the hnear dimensions of C Then the PfafRan 
is the square root of the determinant of a 2dN x 2dN antisymmetric matrix A in the form 
of a direct product of 2d x 2d matrices indexed by 



A(ni, ...,nd\ni, ...,nd) 



^ 



-W 



A(ni, Ui, nd\ni, Ui + 1, na) = 







A{ni, Ui, nd\ni, Ui - 1, Ud) 



^ o^ 



V 



A(mi, ...,md\ni, ...,nd) 



-U, 
otherwise 



(7) 



where is the d-dimensional zero matrix, W a d x matrix with elements — ooij+dj 
and \Ji a. d X d matrix whose only nonzero element is {Ui}-^ = 1. 

With periodic boundary conditions, matrices in (7) are invariant when indices m^, rii are 
changed by Ni. Then the matrix is block-diagonal in the Fourier space. This leads to the 
following expression for the partition function 



Z= VdetA 



Ni 



n n 

ni=l nd=l 



det 



V 



B 
-B* 



1/2 



= n n det|Bi, 

ni=l nd=l 

where B = B(27rni/A'"i, 2'nnd/Nd) is a d x d matrix with elements 



Bij = 0Ji,j+d, i 7^ j 

= + e2-WAf^ i = j, i,j = l,...,d, 



(8) 



(9) 



and B* is the complex conjugate of B. In this way, we obtain the per-site free energy 
fd= hm N'HnZ 

= ^ d<^i--- de,in\D,{e,,---,e,)\, (lo) 

where 

Dd{eu---,ed) = det{B{eu---,ed)}, (ii) 

Particularly, for d = 2 and 3, we have 

D2i9i, 62) = ^1234 + ^24e*'^ + ^ise^'^ + e*(^^+^^\ 

^3(^1, O2, 63) = ^123456 + tU2356e'^' + ^^13466*^' + ^^12456^^' 

Here, the w's with 4 or more indices are those defined by (2) - (4b). Particularly, we have, 
in addition to those given in (4a) and (4b), u;2356 = ^25^^36 — 1^26^^355^^1346 = i^i4'^36 — i^i6'^34- 
Explicitly, Dd{9i, ■ ■ ■ , 6'^) is a linear combination of terms in the form of Cjj...fce**^^*^^^" ~'"^'=\ 
where the coefficient Cij...^ is the determinant of the matrix W with the i,j, kth rows and 
columns deleted. By taking cuij = ^/ZiZ], Zi = Zi+d for which the free-fermion weights with 
4 or more indices vanish identically, these results reduce to those of [^. 

V. THE CRITICAL BEHAVIOR 

For d = 2 the vertices are those of a six-vertex model shown in Fig. 4 with the weights 

fil = 1, Q2 = 1^1234 = 1^3^24 — 

Q3 = UJi3, Q4 = LJ24 (13) 

^5 = 1-^14) ^6 = '-^23- 

Note that the free-fermion condition is now Q1Q2 = ^3^4, — and differs slightly (cor- 

responding to the negation of ^2) from the usual form ||^ Q1Q2 = —^^3^24 + ^^5^26 for which 
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the partition function is defined by (5) without the loop factor (—1)^. Rewriting (12) and 
carrying out one integration, say, over 62, we have 

/2= r dOi r In 1^2 + ^ze''' + ^^e''^ + n^e'^''+'^^\ 

[Zn)'^ JO Jo 

= — r^rf^ilnmaxjlfia + ^^se^^^l, 1^4 + ^^ic^^^ ||. (14) 
27r JO I J 

This free energy reduces to that of usual free-fermion 6- vertex model after negating Q2 [H- 
The equivalence of the two free energies (with and without the loop factor) is a property 
unique to two dimensions. 

When one of the two absolute values in (14) is larger than the other for all 61, then one 
can carry out the remaining integration and obtains 

/2 = max{|lli|, 1^2!, I^^sl, 1^4!}, (15) 

so that the system is in a "frozen" state. Now the two absolute values in (14) are linear 
in cos6'i. As a consequence, when perturbed from the frozen states (15), the two absolute 
values will cross only at ^1 = or vr, near which one of the two absolute values is always 
larger. This leads to the two critical conditions 

1^2 i^al = |^^4 if^il, (16) 

at which /2 becomes singular. Near the critical point say, Q2 + ^3 = + ^1 for all Qi > 0, 
for instance, let 

t = Q2 + ^3 — — ^ small positive. (17) 



Then, the singular part of /2 is given by 



{Rising = TT n , n ie. ' i^^) 

271 J -a t) iZd + iZl 6**^1 



2tT J-a{t) 

where a(t) > is determined from 



^4 + f^ie^^i 



1^2 + ^^se^^l = |1^4 + ^^ie^°|. (19) 



It is easy to verify that a. ~ and for small t the integrand in (18) is of the order oit^Q\. 
It follow that 

{/2}sing ~ t'-l\ (20) 

The same critical behavior is deduced if any of the weights, say, Q2, is negative. 

The example old = 2 serves to illustrate the origin of the critical behavior and the critical 
point. After carrying out one-fold integration, ji is given in the form of an integration of 
a competition of the logarithms of two absolute values as in the second hne of (14). Then, 
the critical point is the point at which the two absolute values are equal and the integrand 
in (14) switches from one absolute value to the other for small deviations from the critical 
point. This creates a singularity in the free energy. We shall call this the switching property 
of the integrand. A prerequisite for switching to occur is that the ratio of the two absolute 
values is an extremum in Q\ at the critical point. For d=2 this can occur only at 6'i = or 
TT. This leads to a singular part of the free energy in the form of (18), where the integration 
is taken over a region whose boundary is determined by setting the two absolute values 
equal at fixed small deviation t from criticality. It follows that the singular behavior is that 
of this integral, namely, tV{t\ where Vit) ~ ^ (for d = 2) is the volume of 7?.. 

Applying the same argument to c? > 3, one first carries out one-fold integration over, say, 
, yielding in the form of an integration of a competition of the logarithms of the absolute 
values of Aj{Qa)^ j, the collection of terms in hnear in e*^^ , and Bj{6a), the collection 
of terms in Da independent of e*^^ . Then can be singular at points in the parameter space 
at which the two absolute values become equal and possess the aforementioned switching 
and extremum property. Namely, there exist angles 9ao, Oi ^ j, such that 

1^,(^,0)1 = |5,(^ao)|, t^O (21) 

and 

\Aj{e^)\ > \B^{e^)i t<o 

\Aj{ea)\ < \Bj{e^)\, t^o+ (22) 
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for 9a in a neighborhood of Oao-, where t is the deviation from the critical point. (Here, the 
roles of Aj and Bj can be interchanged.) Since the integrand of the free energy switches 
from \Aj\ to \Bj\ for t = 0+, a singular part in the form of (18) appears in the free energy, 
which is the integration of In \Aj/Bj\ over a region TZ bounded by 

\Aj{9a)\ = \Bj{9a)\, t = small fixed. (23) 

The critical condition is now given by (21) and the critical behavior is tV{t), where V{t) is 
the volume of the {d — 1) -dimensional region TZ. 

To determine V{t), we expand In \Aj/Bj\ about 6*00 and t = 0, obtaining generally the 
leading behavior 

ct + F{Aea), (24) 

where c is a constant and A9a = Oa — Oao- Here, due to the extremum property, F is a 
quadratic form in variables A6a, which is either positive or negative definite. This implies 
that the region TZ has linear dimensions A6a ~ t^^"^, and hence V{t) ~ t('^^i)/2. The singular 
part of the free energy behaves as 

{/4sing~t('^+')/'. (25) 

This is the critical behavior found in the model = [§. The meaning of the singular 
behavior (25) for d = odd is that the [{d + l)/2]-th derivative of fd is discontinuous. 

However, the singular behavior (25) may not be the dominate one. If the {d — 1)- 
dimensional integral /ln|y4j| dominant in (21) - (23) is also singular at t = 0, then the 
critical behavior of the free energy at t = is that of a (rf — l)-dimensional model. By 
the same token, the critical behavior can be reduced further to that of lower dimensions. 
Examples are given in the next section. 

VI. EXAMPLES 
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A. Example 1 



If lines in the vertex model are nonintersecting, then we have the model considered in 
§] with 

/^d(^l,■■■,^d) = l + E^.e*'^ (26) 
1=1 

There is no loss of generality in supposing \zi\ < 1. Carry out the integration in (10) over 
6i, say, after rewriting Dd as e~*^^ + zi + ^2^*^^ + ■ ■ ■ + zae'^^'', where we have renamed 6a — di 
as 6a- If any of the Za is negative, one replaces 6'^ by vr + 6a and obtains Ai = 1 and 
Bi = \zi\ + 1 2^2 1 e*^^ + ■ ■ ■ + \zd\e^^''- It is then straightforward to verify that the extremum 
property occurs only at 6ao = 0, a = 2, ...,d. The critical point occurs at 

\zi\ + \Z2\ + ■ ■ ■ + \zd\ = I, (27) 

and the critical behavior of the free energy is 

^(d+l)/2_ 

B. Example 2 

. Consider the 3-dimensional model with weights 

1^36 = Z, UJ24 = Zi/ Z 

i^ie = = z 
0^34 = 1^35 = 1 

0^14 = = a;i5 = (28) 
for z and Zi real, implying, using (4a) and (4b), 

^^123456 = Zii (^1245 = 

'^1346 = 1^2356 = —Z. (29) 

By appropriately changing integration variables, the per-site free energy (10) assumes the 
form 
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1 rln r27r pItt 

h = . / dOi / (192 / dOz In bi + ze'^' + ze'^^ + ze'^^'^^^^ + e^(^^+^2+^3) I . (30) 

(27r) ^0 ^0 ^0 

After carrying out the integration over ^3, one obtains — 1, — zi -\- ze^^^ + ze^^^ + 
with 

l^sl^ = {zx + 'izf + 2z{z + ^i)(cos^i + cos 6*2 — 2) 

+ 2zzi[cos(^i + Q2) - 1] + 2^2[cos(^i - Q2) - 1]. (31) 

It can be verified that \B-^ possesses one extremum at ^10 = ^20 = in the regimes Zi/z > 
— 1/3 and zi/ z < —3, which, by setting the extremum equal to one, leads to the critical 
point 

\zi + 3z\ = 1. (32) 
Define t = \zi + 3z\ — 1. Then one finds 



for zi/z > —1/3, and 



/3 = 0, t<0 

{/ajsing ~ t', t = 0+ 

/3 = ln|^i|, t>0 

{/alsing ^t^, t^O-. 



(33) 



(34) 



for zi/z < —1/3. In addition, another extremum occurs at cos6'io = cos6'2o = —{^i + z)/2zi 
in the regime |1 — Szi/z\ > 2, which leads to the critical point 

{zi - zf = zi. (35) 

Define t — \zi — zf — \zi\. Then one finds the behavior (33) for zi/z < —1/3 and (34) for 
zi/z > —1/3. These results lead to the phase diagram shown in Fig. 5. 



C. Example 3 

Consider another 3-dimensional model with weights 
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^15 — UJ^Q — —1 

^lA = a;25 = 

<^24 = <^26 = <^34 = <^35 = 0^16 = 1 (36) 



implying, 



'^2356 = '^1346 = ~1 
<^1245 = — ^ 

(^123456 = 2z + l. (37) 



The per-site free energy is 

\ i-2n r-2n i-2-K 



fZn fZn rziz 

(27r) -^0 Jo Jo 



(38) 



After carrying out the integration over ^3, one obtains A3 = z — e'^^^^~^^^\ = 2z + 1 — 
^1(61+62) _ giSi _ ^102^ leading to the critical point 

z^l, (39) 

at 6*10 = 6*20 = 0. Define t = z — 1. It can be verify that lA^] < IB^] for t > 0, and a switching 
from I S3 1 to \As\ occurs in a small regime in the neighborhood of ^10 = —^20, t — 0—. Hence 
the free energy contains a singular part of the order of 0{t'^) for t — 0—. However, it turns 
out that the dominate integral 

7A2 rde^ r d92ln\Bsi9^,92)\, (40) 
(27r) Jo Jo 

is also singular at t = 0. As in (14), further integrations of (40) lead to the critical behavior 

/3=ln|2z + l|, t>0 
= ln|2^ + l| + ^|t|3/2, t^O-. (41) 

This critical behavior is verified by numerical integrations whose results are shown in Fig. 
6. 
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VII. SUMMARY 



We have solved exactly a free-fermion vertex model in d dimensions. The vertex model is 
described by lines running in a preferred direction with the weights of vertices of intersecting 
lines prescribed by the free-fermion condition (2). The partition function (5) contains a 
fugacity —1 for each loop of lines and therefore consists of positive and negative weights. 
The per-site free energy is evaluated exactly in (10). The critical point is found to be given 
by (21) and the critical behavior given by In exceptional cases this critical behavior 

is modified with d replaced by an integer lower than the actual dimensionality. 

This work has been supported in part by NSF Grants DMR-9313648 and INT-9207261. 
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FIGURES 

FIG. 1. Decomposition of a vertex configuration of 2r incident bonds into r! line configurations. 

FIG. 2. A city of 2d points for d = 3. 
FIG. 3. The correspondence between vertex configurations on £. and dimer coverings on 

FIG. 4. Vertex configurations and weights for a 6-vertex model. 

FIG. 5. Phase diagram of the three-dimensional model (30). Regimes I and II are in frozen 
states with /s = and fs = ln\zi\, respectively. 

FIG. 6. Results of numerical integrations of dfz{z)ldz for the three-dimensional model (38). 
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